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Laboratory experiments, large-scale computer simulations and observational cosmology have be¬ 
gun to make progress in the campaign to identify the particle responsible for gravitationally-inferred 
dark matter. In this contribution we discuss the dark matter density profiles in the cores of nearby 
galaxy clusters and estimate the gamma-ray flux expected for MSSM dark matter over a range of 
neutralino masses. 


I. INTRODUCTION 

Galaxy clusters are dark matter-dominated objects. 
Since they are believed to constitute a nearly fair sam¬ 
ple of the matter content of the universe [lj, WMAP 
results suggest that they are comprised of roughly 
15% baryons and 85% non-baryonic dark matter j3|- 
Numerical experiments which simulate their forma¬ 
tion via the hierarchical assembly of cold dark matter 
(CDM) halos suggest that their density profiles are 
adequately described by a pair of power laws and a 
transition radius 0 a 1 . ' Indeed, the dark matter 
density profile in the centers of relaxed clusters shows 
power-law behavior to scales as small as ~ 10 kpc 
|6, 3, |8!j. If the density profile remains a power law 
to small enough radius, the central density may be¬ 
come large enough that dark matter self-annihilation 
produces a gamma-ray flux observable with current 
instrumentation. 

The neutralino is perhaps the leading candidate for 
the dark matter particle {9|. The observational sig¬ 
nal for neutralino annihilations can be quite spectac¬ 
ular. Two of the annihilation channels result in the 
production of monochromatic gamma rays, unlikely to 
be confused with other astrophysical processes. Other 
channels can lead to substantially more gamma rays, 
although since they produce a continuous spectrum 
they are more difficult to distinguish from high en¬ 
ergy processes such as shock heating. 

Most searches for the gamma ray signature of dark 
matter annihilation have centered on local dark mat¬ 
ter concentrations such as the Galactic center [T(j and 
halo 00 and Local Group dwarf spheroidals 0 - 
However, if the centers of galaxy clusters are cuspy 
to sufficienty small scale, they too may be observable 
with the current generation of gamma ray telescopes. 
In this contribution we calculate the expected annihi¬ 
lation signal from clusters at distances less than ~ 100 
Mpc for core density profiles determined through X- 
ray observations. 


ton, Z°, and neutral Higgs bosons, and that it pro¬ 
vides the bulk of the dark matter in the universe 
0 0 The neutralino, a Majorana fermion, self- 
annihilates in the early universe until the annihilation 
rate is exceeded by the Hubble rate. The neutralino 
relic density Ul x depends upon its exact composition 
and the presence of resonances, but it can be deter¬ 
mined approximately using 0 


Ul x h 2 


3 x 10 27 cm 3 s 1 
(av) 


(1) 


where h is the Hubble parameter and (av) is the ther¬ 
mally averaged annihilation cross section. Assuming 
neutralinos are by far the cosmologically dominant 
dark matter species, then ~ (1 — fb) H m , where fb 
is the matter baryon fraction and Ul m is the total mat¬ 
ter density. For h = 0.67 and fb = 0.15, the total neu¬ 
tralino annihilation cross section is (av) = 3 x 10 -26 
cm° s . 

Neutralinos can annihilate directly into a pair of 
monochromatic gamma rays [77], a gamma and a 
neutral Z boson [7Z 0 ], or into a gamma continuum 
through a plethora of hadronization processes [7(h)]. 
The dominant “channel” is hadronization, with a re¬ 
sulting gamma ray spectrum conveniently approxi¬ 
mated by 0: 


d 2k = j^r dx a-*) 2 (2) 

dE 1 4m x J Xmin x 3 / 2 (x 2 — rj 2 ) 1 ! 2 {) 

Here x = E v /m. x , rj = m v /m x , and x m i n = E^ + 
m x rj 2 /AE X . The neutral pion and neutralino masses 
are m n and m x . The 77 and 7Z 0 spectra are given by 


dN i = 2 _ s ( 1 _ m x \ 

dE 7 E \ E 7 ) 

and 


( 3 ) 


II. NEUTRALINO ANNIHILATION 

We will assume throughout this contribution that 
the lightest supersymmetric particle is the neutralino, 
a linear superposition of the superpartners of the pho- 


dNj_ = J_ s f 1 _ m x( 1 - (m z o/2m x ) 2 ) 
dE 1 E \ E 1 

respectively 000 . Estimates for the cross sec¬ 
tion of these last two channels vary over several or¬ 
ders of magnitude. We use the following values for 
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FIG. 1: Gamma ray annihilation spectra for four different 
neutralino masses. The spectra shown are, left to right, 
for 10, 30, 100 and 300 GeV/c 2 . The inset shows the indi¬ 
vidual contributions from the 77 , 7 Z 0 and 7 (h) processes. 


our analysis: (av) 7 ( h ) = 3 x 10 -26 cm 3 s _1 , (cru) 77 = 
(crv)j’z 0 = 0.01(cru) 7 (h). These spectra are shown in 
Figure 0 for four values of the neutralino mass. 


III. GAMMA RAY LUMINOSITY OF 
CLUSTERS 

The annihilation radiation luminosity of a galaxy 
cluster is given by 


Nj ( av) 
47TTO^ 


p 2 dV 


( 5 ) 


where p is the matter density of neutralinos, and the 
integral is taken over the entire cluster volume. For 
a dark matter core profile following p = kr a , this in¬ 
tegral diverges toward the center for a < —1.5. Al¬ 
though astrophysical processes will probably impose a 
density cut-off, significant gamma ray signatures may 
still result. 

Figure [3 shows the inner logarithmic density slope 
for a sample of seven relaxed clusters observed with 
the Chandra X-ray Observatory || 23. The slopes 
generally span the range — 1 < a < —2, suggesting 
that gamma ray luminosities may indeed be quite high 
at the center of some clusters. 

The largest astrophysically possible neutralino den¬ 
sity is determined by the free-fall timescale at the cen¬ 
ter, with the density determined by equating the free- 
fall rate 


( 6 ) 
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FIG. 2: Inner logarithmic density slopes for seven relaxed 
galaxy clusters 6, Q. The gray band, bracketed by the 
NFW and Moore profiles, represents the prediction of stan¬ 
dard CDM cosmology. 


to the annihilation rate 


ta 1 =n x(° v ) ( 7 ) 

since neutralinos cannot annihilate faster than they 
are supplied to the high density regions. This implies 
a free-fall-limited maximum density 


PmaxJ — 1-2 X 10 ( * 2 ^ M© pc (8) 

ylOGeV/c ) 

However, dynamical heating of neutralinos by the for¬ 
mation of supermassive black hole binaries may re¬ 
duce this value significantly, to about 10 5 M 0 pc - 3 
19!]. (Other processes such as gravitational heating 
pOj may be less restrictive.) Figure 0 illustrates the 
peak central density required for a galaxy cluster to be 
rendered detectable by GLAST and HESS as a func¬ 
tion of cluster distance. The cluster central density 
profile used is a = —1.5 to a cut-off radius, where it 
remains at the limiting density. The horizontal black 
line indicates a crude estimate of the upper limit of 
the density imposed by the formation of supermassive 
black holes during the growth of the central cluster 
galaxy G3- This figure shows that if central densities 
reach high, but plausible values, then annihilation ra¬ 
diation from neutralino dark matter may be detectable 
to distances of tens to hundreds of Mpc. This radi¬ 
ation would carry important information about the 
very centers of these objects. 


IV. A TALE OF TWO CLUSTERS 

It makes sense, in light of these results, to ask the 
question: Could any real clusters be detected in their 
annihilation radiation? To address this question we 
use two examples. The first is the Virgo Cluster at 
a distance of 17 Mpc, with a central density slope of 
— 1.3 j^J. For the second example we use an Abell 


tj 1 = 32Gp/3n 
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FIG. 3: Dark matter central density required for detectability as a function of cluster distance. The upper three diagonal 
lines represent GLAST sensitivity limits for m x = 300, 100, 30 GeV, the lower two for HESS sensitivity limits at 300 
and 100 GeV. 


2029-like cluster (a = —2; 0) at a distance of 100 
Mpc. The annihilation radiation flux from a cluster, 
in photons s -1 cm -2 sr -1 , is 


neutralino mass is low; HESS may be able to detect 
it for a wider range of distance and particle mass. 

V. CONCLUSION 


F 1 = 2.9 x 10 


-7 


(<TV) 


■( 


3 x 10 -26 cm 3 s _1 
\ -2 


m y 


1 


dfl / AT 


yiOGeV/c / ^ Jn je 

[ ( p y di 

Jlos l 10~ 24 g/cm 3 J lkpc 


( 9 ) 


where is the beam size. Figure0]shows the expected 
annihilation radiation flux from each object as seen by 
GLAST (for 30 < m x < 300 GeV/c 2 ) and HESS (for 
100 < m x < 300). The GLAST and HESS detection 
thresholds are shown for comparison, as is the ten¬ 
tative HESS detection of Sgr A* !22]. In both cases 
we have adopted p m ax = 10 5 M 0 pc -3 . The fluxes 
have been calculated from the surface brightness inte¬ 
gral over a beam size of 10 -5 sr. GLAST may detect 
an annihilation signal from a steep-core cluster if the 


We have shown that, for plausible values of the cen¬ 
tral dark matter slope and for sensible astrophysical 
constraints on the peak density, neutralino annihila¬ 
tion radiation from the centers of galaxy clusters may 
be detectable using current and near-future gamma 
ray telescopes. These signals would not only provide 
spectacular, if indirect, evidence for the existence of 
dark matter, but they would also carry information 
about conditions at the centers of galaxy clusters. 
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FIG. 4: Expected gamma ray fluxes from Virgo (left at 17 Mpc) and a more distant, but steeper-profiled, cluster at 100 
Mpc. For reference the GLAST and HESS detection thresholds are shown, as is the recent detection of Sgr A* [22]. 
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